To what extent does the self-consistent mean-field exist? 
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A non-convergent difficulty near level-repulsive region is discussed within the self-consistent mean- 
field theory. It is shown by numerical and analytic studies that the mean-field is not realized in the 
many-fermion system when quantum fluctuations coming from two-body residual interaction and 
quadrupole deformation are larger than an energy difference between two avoided crossing orbits. 
An analytic condition indicating a limitation of the mean-field concept is derived for the first time. 
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The Hartree-Fock mean-field theory has been greatly 
successful in describing many-fermion systems. Level 
repulsion is a universal phenomenon appearing in the 
mean-field description of molecular, atomic, biological 
and nuclear systems [H, 0, S 0| • In the development of 
nuclear structure physics, there have been many discus- 
sions on the applicability of cranked mean-field theory 
near level-repulsive region [EIQ. An argument on how 
to remove certain spurious interaction and to construct 
diabatic orbits seems to be reasonable in cranked mean- 
field [7j, |8j, because two crossing orbits interact not at 
a given angular momentum but at a given angular fre- 
quency. However the above argument may not be sim- 
ply extended to the deformation constrained mean-field 
theory and the interaction between two potential energy 
surfaces (PESs) might not be necessarily regarded to be 
spurious since the Hamiltonian and deformation opera- 
tors do not commute with each other. Although the non- 
convergent difficulty in the mean-field calculation was re- 
alized numerically [7j], to our knowledge, there has been 
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no discussion by deriving analytic expressions to clearly 
demonstrate why and how the competition between one- 
body potential and two-body residual interaction plays a 
decisive role in breaking the concept of mean-field near 
level-repulsive region. For this purpose, this Letter will 
give an analytic condition to indicate the applicability of 
mean field near level replusion. 

The quadratic constrained Hartree-Fock (CHF) equa- 
tion with Gogny D1S interaction @ has been solve d by 
using three dimensional harmonic oscillator basis [lfj |. 
The triaxial oscillator parameters are determined by min- 
imizing the energy of HF solution. In the calculation 
of PES, the same oscillator parameters as those of the 
ground state are used to trace an evolution of the ground 
state configuration, and to make the single-particle (s.p.) 
level-crossing dynamics transparent. Figure [Ha) shows 
that there appears a missing region in the PES of nucleus 
66 Se. At an edge of the PES ((Q) = q « 178 fm 2 ), the 
CHF calculation meets a difficulty of non-convergence 
even one decreases the deformation quite a little, and no 
matter how much effort one makes to get convergence. 



FIG. 1: Ground state PES of nucleus 66 Se (a); diagonal com- 
ponents of proton (b) and neutron (c) densities as a function 
of iterative number for a non-convergent case; neutron s.p. 
energies near Fermi surface (d). 



After the missing region, a continuously-connected PES 
passing through another local minimum (q w 100 fm 2 ) is 
obtained. In order to understand what prevents the CHF 
calculation from convergence, the diagonal components 
of proton (b) and neutron (c) densities as a function of 
iterative number are depicted for a non-convergent case 
with fj, = 150 fm 2 . Here, /i is a control parameter of 
quadratic constraint given by —w((i— (Q))Q, and is used 
to assign deformation in the missing region instead of q. 
A representation where one-body CHF hamiltonian h^ n ' 
is diagonal is used with n being iterative number. From 
FigOJb), one ma y observe that the expectation values of 
proton density converge to for unoccupied orbits and 
1 for occupied orbits. For the case of neutron density 
(c) , there appears a similar situation for most single- hole 
and single-particle states, except for two specific orbits 
labeled as (— , — )a and (— , — )a+i, where (— , — ) denotes 
parity and signature quantum numbers and A the num- 
ber of occupied orbits. The two specific orbits are inter- 
acting and lying just below and above the Fermi surface 
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as shown in Fig.QJd). One may expect that the two spe- 
cific orbits play a dominant role in preventing the CHF 
calculation from convergence. 

In what follows, we will focus on how the non- 
convergent difficulty appears as a result of the micro- 
scopic dynamics. Suppose there exists a convergent CHF 
state at q = qo. We will discuss whether the CHF calcu- 
lation with slightly different deformation at q = qo — Aq 
goes to convergence or not. When one takes the CHF 
state at q = q as a starting trial function to proceed an 
iteration at q = q — Aq, the matrix elements of the CHF 
Hamiltonian in the first iteration between two specific or- 
bits in a qo-representation of using the self-consistent s.p. 
states {ek{qo), <fik(qo)} at qo are written as (for detailed 
analytic derivation, see Ref. (Tlj |) 
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where a, (3, ■ ■ ■ refer to a set of s.p. states used for numer- 
ical basis. In the go-representation, the matrix elements 
of deformation operator in Eq.([T]) are expressed as 
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The 2x2 truncated Hamiltonian of h^(qo 
first iteration is expressed as 
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Making the following analytic discussion simple, we ex- 
ploit such an approximate expression that a contribution 
to the mean-field from the (A — l) number of hole-states 
is independent of small deformation change A/x, and of 
number of iterations, which is well realized in our numer- 
ical calculation discussed above. That is 
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Under the above assumption, one may explore the non- 
convergent dynamics governing the CHF calculation in 
terms of the 2x2 truncated CHF Hamiltonian. 

Although a set of s.p. eigenstates {e^ (go — 
Aq), (p^\qo — Aq)} is numerically obtained by diagonal- 
izing the full CHF Hamiltonian (qo — Aq) , a charac- 
teristic feature of interacting orbits at the nth iteration 
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FIG. 2: For a given non-convergent case the off-diagonal 
Hamiltonian and difference of diagonal Hamiltonians (a) and 
the value of a (n) & (n) (b) as a function of iterative number. 



may be understood by the truncated Hamiltonian and 
2x2 unitary matrix 



a {n) h (n) \ / 

dto C W J I y A +x(qo) 




= u ir 



<PA(qo) 
¥>A+i(qo) 



(5) 



In each diagonalization, A and ^4 + 1 are always used to 
label the s.p. states in an energy increasing order. To 
show a decisive effect of relative phase a^b^ on the 
non-convergence of the CHF calculation, we introduce 
four inter-dependent parameters in rather than a 

single independent parameter. The mixing parameters 
and measure a degree of interaction between 
two specific orbits. Here and hereafter, ip^ (qo — Aq) and 
h^ n \qo — Aq) at the nth iteration are simply expressed 

as ip^ an d h^ n \ 

In the go-representation, one gets the relation between 
two successive CHF Hamiltonians as 
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and the anti-symmetried two-body interaction 
va+iaaa+i is defined in the go-representation as 

VA+1AAA+1 



Va+Ic (<1o)<Pa 7 {qo)v ai 06<Pi3A (q )(fsA+i (<?o) ■ 



(8) 



Similarly, the (n + 2)th off-diagonal CHF Hamiltonian is 
related with the nth through 

*<%a+i = h AA+i + {a (n+1) 6 (n+1) - afr-W-V} 



VA 



-iaaa+1 - {AA<" +1 > - AA^}Q 



(9) 



Numerical values of h 



in) 
a,a 



^a+i,a+i an( i h A ? A+1 are 



shown in Fig. [DJa) for the non-convergent case with 
(j, = 150 fm 2 . Here the matrix elements are calculated 
in the go-representation at the edge point of convergence 
go = 178 fm 2 . From this figure, one may observe that 
both components exhibit a staggering property around 
some averaged values. It should be noticed that the off- 
diagonal hamiltonian h AA , j changes its sign from iter- 
ation to iteration. According to the perturbation theory 
for 2x2 Hamiltonian, a sign of a^b^ of the lower 

' l A+l,A+ll- 111 



state is given by that of (h^ A \ +1 ) / (h 

Fig-El ), a numerical value of a^ n 'b^ n ' is shown as a func- 
tion of iterative number. Comparing Fig. [21 a) with (b), 
one may recognize that the above discussion of perturba- 
tion theory is well justified by the numerical calculation. 
Since the sign of a changes by iteration and iter- 

ation, one may deduce such a conclusion that the prop- 
erties of the upper state (— , —)a+i and the lower state 
(—,—)a ar e inter-changed from one iteration to the next. 

As seen from Eq. ©, h AA +i becomes large because 
of the different signs between a^ n 'b^ n ' and a*" -1 ^" -1 -*, 
whereas ft^^+i becomes small at the (n + 2)th iteration 

because of the same sign of a(" +1 )&(™ +1 ) and a^™" 1 ^™ -1 ) 
in Eq. Since the quantum fluctuations coming from 
two-body residual interaction and quadrupole deforma- 
tion become small in one iteration and large in the next, 
there appears the staggering property and they could not 
be approximated successfully by the one-body mean-field 
potential. Physically one could understand the above sit- 
uation as follows: two mean-fields, one characterized with 
tp A (go) and the other with tpA+i(qa), interact too strongly 
by the two-body residual interaction to be approximated 
by a single mean- field. Although the above phenomenon 
in Fig. [5] has been known as "ping-pong" [7j, it should 
be noted that the present analytic understanding is given 
for the first time. 

The above analytic and numerical results clearly in- 
dicate that the sign of a (")&(") between two successive 
iterations is important whether or not the self-consistent 
CHF theory would be applicable. Our next task is to 
make clear an analytic condition, which could tell us 
an existence of the mean-field. The preceding discussion 
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FIG. 3: The difference of s.p. energies and quantum fluctua- 
tion. 



clearly indicates that the concept of CHF mean-field is 
applied successfully, provided there holds a condition 
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whose most important part is evaluated for a case with 
n = 1. With the aid of Eqs.((6|) and ([9j) , it is easily shown 
that Eq. (fTU|) is satisfied for any n provided there holds 
Eq. (fT0| with n = 1. After some calculation, one may 
derive a relation 

r 1.(2) i 



sign -j ^ = sign{eA+i(9o) - e>l(?o) ~ v a +iaaa+i 

I h'A.A+l 



-2wQ 2 A<A+1 - wAfi{Q AtA - Q a +i,a+i) - 0(b)}. (11) 

Here 0(b) contains a small mixing parameter fo( 1 ), and is 
expressed as 

0(b) = -2b (1)2 (va+iaaa+i + 2wQ\ A+1 ) + 
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Condition (|10p for n = 1 is then expressed as 

£A+l(go) - £A(go) > VA+1AAA+1 + 2wQ AA+1 

+wA/j,(Q a ,a - Qa+i,a+i) + 0(b). 



(13) 



Here the r.h.s. of condition ([TBI with neglecting the 
higher-order effect 0(b) is the quantum fluctuation 
expressed as va+iaaa+i + 2wQ A A+1 + wA^l(Q AiA — 
Q a+i,a+i)- The physical meaning of condition (|13| is 
clear: the two-body correlation between nucleons can be 
successfully incorporated into the mean field when the 
energy difference between two interacting orbits is not 
smaller than the quantum fluctuation coming from two- 
body residual interaction and quadrupole deformation. 
Similar derivation gives the opposite condition 

e A+i(go) - £A(go) < v A+ i AAA+ i + 2wQ 2 A A+1 
+wAfi(Q A ,A - Qa+i,a+i) + 0(b), (14) 
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FIG. 4: The ratio of off-diagonal Hamiltonians between the 
first and second iterations for each given /x. 

which states the breakdown of mean field. 

Figure [3] shows the l.h.s of condition (fl~3"|) as well as 
the r.h.s. with neglecting the higher-order effect 0(b). It 
is clear that our numerical results show the feasibility of 
the analytic condition (fl3|) . Namely, both sides of condi- 
tion (|13p take almost the same value, when q reaches to 
the edge points of the convergence. Figure [4] depicts the 
ratio h^\ + i/h^' A , v the l.h.s of Eq. (jTUJ) with n = 1, 
as a function of quadrupole parameter /i. One may see 
that the ratio is kept positive in convergent region, and 
reaches zero near the edge of it. When one slightly de- 
creases \x from the edge point, the ratio becomes negative 
where the convergence of CHF calculation is not achieved 
and condition (fl~4f is satisfied. 

The present work does not include the pairing cor- 
relation. Including the pairing correlation might make 
the situation so far discussed more complicated because 
of many dynamical competitions that exist not only be- 
tween the p/i-type two-body residual interaction and the 
HF potential, but also the pp-type two-body residual in- 



teraction and the pairing potential; their cross effects also 
must be considered. Moreover, the two mean fields char- 
acterized by different configurations are mixed up by the 
pairing correlation, and the uv factor introduced by the 
BCS theory obscures the concept of the configuration. 
Some numerical evidence of the above non-convergent 
feature in CHF-Bogoliubov (CHFB) theory has been dis- 
cussed in Ref. flo| and its analytic derivation on an ap- 
plicability of the CHFB theory near level crossing region 
is on progress. 

Summarizing the present Letter, we have shown that 
the concept of CHF mean-field breaks down near the 
level-repulsive region when the condition (|13[) is not 
satisfied. Deriving the analytic condition, we make it 
clear that the competition between one-body potential 
and quantum fluctuations mainly coming from two-body 
residual interaction plays an important role whether the 
self-consistent CHF mean field is realized or not. Since 
the above result is obtained by using the constrained op- 
erator Q put by hand, further studies are needed whether 
an introduction of "diabatic" orbits @, 0] really remedies 
the present difficulty or not. That is, the system does 
not like to be elongated nor contracted along a given di- 
rection of quadrupole deformation any more, but it likes 
to develop toward a direction chosen by itself. Further 
microscopic investigation is needed to answer the very in- 
teresting conclusion in Ref. 0] by introducing dynamical 
constrained operators based on the self-consistent collec- 
tive coordinate (SCC) method [H,[l3j], because a diabolic 
point related to the level crossing indicates an existence 
of missing degree of freedom. 
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